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Abstract 

There are good reasons why neutrinos could be Majorana particles, but there exist 
also a number of very good reasons why neutrinos could have Dirac masses. The 
latter option deserves more attention and we derive therefore analytic expressions 
describing the renormalization group evolution of mixing angles and of the CP phase 
for Dirac neutrinos. Radiative corrections to leptonic mixings are in this case en- 
hanced compared to the quark mixings because the hierarchy of neutrino masses is 
milder and because the mixing angles are larger. The renormalization group effects 
are compared to the precision of current and future neutrino experiments. We find 
that, in the MSSM framework, radiative corrections of the mixing angles are for large 
tan/3 comparable to the precision of future experiments. 
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1 Introduction 



One of the most important open questions of neutrino pliysics is whether neutrinos are 
Dirac or Majorana particles. Prom a theoretical perspective, large Majorana mass terms 
appear quite naturally for the right-handed neutrinos, since they are complete gauge sin- 
glets. This leads directly to the best investigated (and therefore most popular) explanations 
for the huge ratio of observed mass scales in the see-saw mechanism [1-6] . In its simplest 
form, neutrino masses get suppressed by a factor vew/M^, with t;Ew denoting the vacuum 
expectation value (VEV) of the Higgs boson and being the scale at which B — L sym- 
metry (baryon minus lepton number) is assumed to be broken. However, it is important 
to keep in mind that the suppression factor v/M^ (with M^, now being the GUT scale or a 
related scale) arises rather generally whenever neutrino masses arise from integrating out 
heavy degrees of freedom with mass M*. This statement holds independently of the nature 
of neutrino masses, in particular both for Dirac and Majorana masses. Indeed, there exist 
a couple of appealing models where small Dirac masses are explained in this way by using 
extra, heavy degrees of freedom, or by relating the Yukawa coupling to the ratio of 
gravitino mass (or other soft masses) and GUT (or compactification) scale [7-12]. An- 
other possibility is using localization in extra dimensions, and explaining the suppression 
by a small overlap of the corresponding zero-mode profiles along extra dimensions (see, 
e.g., [13-15]). Further support for Dirac neutrinos was found in certain orbifold compacti- 
fications of the heterotic string, where it is difficult to obtain the standard see-saw [16]. For 
recent overviews and further references of various possibilities of explaining small Dirac 
masses see [17-19]. 

Cosmological arguments do not give a preference for Dirac or Majorana masses either. 
For instance, even if one requires the explanation of the observed baryon asymmetry to 
be related to neutrino properties, one finds that successful baryogenesis can work both 
for Majorana [20,21] and Dirac [22,23] neutrinos. Dirac neutrinos evade also constraints 
from primordial nucleosynthesis, since the right-handed degrees of freedom decouple with 
a low temperature so that their energy density is relatively suppressed [24] . The question 
whether neutrinos are Dirac or Majorana particles is therefore one of the main motivations 
for improved neutrino-less double beta decay experiments. Both options should therefore 
be studied seriously until this question is clarified by experiments. 

We investigate in this Letter RG effects under the assumption that neutrinos are Dirac 
particles, and that the small Yukawa couplings are explained by some mechanism operating 
at a high, e.g. GUT or compactification, scale, denoted by Mqut in the following. The 
radiative corrections to the leptonic CP violation has been studied in [25] . We extend this 
analysis to include all leptonic mixing parameters, and derive analytic formulae describing 
the renormalization group (RG) running of the leptonic mixing parameters. The radiative 
corrections are compared to analogous corrections in the quark sector, and we will show 
that generically RG running in the lepton sector is stronger than in the quark sector since 
the coefficients of the renormalization group equations (RGBs) are enhanced due to the 
fact that mass hierarchy is milder and the mixing angles are larger in the lepton sector. 
We compare the size of the radiative corrections to the accuracy of present and future 
neutrino experiments, and find that they are particularly relevant if neutrino masses are 
degenerate and tan /3 is large in the MSSM. 
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2 Analytic formulae 



Using the standard parameterization (see, e.g., [26]) for leptonic (and quark) mixing the 
RGEs for the leptonic mixing angles read 



-C y'^ m\ + 

32 TT^ 777,2 ~ "^1 



ei2 = ^777^ -i f sin(2 ^12) sin^ ^23 + ^(^13) , (1) 



-C y'^ 1 

^13 = 00 ^ 2\ I 2 2T {("^2 - "^1) "^3 cos 5 cos 6*13 sin(2 6'i2) sin(2 6'23) 

+ [7773 — (7772 — "^1) cos(2^^i2) —777^7772] cos^ 6*23 sin(2 6'i3)} , (2) 
h -C yl \m\ -m\ml + {ml -ml)ml cos(2 ^12)] ^ ^,^(n n (n^ 

= -32^ {ml - ml) {ml - ml) sm(2M + m3), (3) 

where the dot indicates the logarithmic derivative w.r.t. the renormalization scale /z, e.g. 
6^12 = d6'i2/dt = /id6'i2/d/i, and 

r 1 , (MSSM) , 

^ ~ \ -3/2 , (SM) . 

Here, we have neglected the tiny electron and muon Yukawa couplings, as well as the 
neutrino Yukawa couplings, against the r Yukawa coupling t/^-. Furthermore, in 6^12 and 
^^23 we only kept the leading order term of an expansion in the reactor mixing angle 9i^. 
We have checked numerically that this is a good approximation for realistic values of 6*13. 

It is instructive to compare Eqs. (^-(jS)) to the corresponding expressions for Majorana 
neutrinos. Technically one obtains Eqs. (0)-® by discarding all terms which depend on 
the Majorana phases in Eqs. (8)-(10) of Ref. [27]. One could thus say that the running 
of the Dirac mixing parameters resembles the running of Majorana mixing parameters 
averaged over the Majorana phases and (p2^- This means in particular that strong 
damping effects for the evolution of 6'i2, as observed for Majorana neutrinos, cannot occur 
in the Dirac case. 

From Eqs. (Q)-©, we can immediately recognize several features of the RG evolution. 
First, for a strong mass hierarchy, the running of the angles is negligible. For 7771 = 0, the 
angles always run less than 1° (except for ^23 which runs more if tan/3 > 40). Second, for 
mi > 0.02 eV, 612 has the strongest RG evolution. Third, as is obvious from Eqs. ((H) and 
(jHl), in the MSSM 9i2 always increases when running downwards while 6'23 increases for a 
normal and decreases for an inverted mass hierarchy. This means that these two angles 



^Stated differently, the running in the Dirac case is roughly half of the maximal running in the Majorana 
case. The factor 1/2 can be understood from the structure of the RGE: in the Dirac case, the mass matrix 
gets rotated by only one term (cf. Eq. ljA.3dp '). 

Am,y = C m,y (YgtVe) + flavor-trivial terms , 

while in the Majorana case there are two terms 



Am^ = C 



7, {y}Y,) ^ {Y^Y,)'^ 



flavor-trivial terms 



with C = -3/2 in SM [28] and two-Higgs models [29], and C = 1 in the MSSM [30,31]. 
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are radiatively enhanced (for normal mass ordering) which may, at least partially, be the 
reason for their large size. Whether ^13 increases or decreases depends on 6. 
The evolution of the Dirac phase 6 is described by^ 

5 = ^(-D^-i + 5(0) + 5(1) + ^(^2^) , (5) 
with the first two coefficients 5^'^'^ given by 

= ^ . ~ f sin(^) sin(2g,,) sin(2g,3) , (6a) 

32 vr^ (777,3 ~ ''^1) ~ ^2) 

= 0. (6b) 
Moreover, the term linear in contains 

= ^,T-2 ~ 2 2Y^ot(Msin(2Msin5+... , (6c) 

which becomes relevant if 61^ is not too small. 

As usual, 5 is undefined for 6'i3 = 0. Clearly, if 5 vanishes for some scale, the (lepton 
sector of the) theory is CP invariant at this scale and thus remains CP invariant for all 
scales. Hence, the statement 5 = cannot depend on the renormalization scale, which 
can also be seen in our formulae. Likewise, if ^13 is zero at some given scale, the theory 
must again be CP invariant for all scales'^. From this we conclude that ^13 can never cross 
zero if we have at some scale ^13 7^ and sin 5 7^ 0. If ^13 approaches zero, then we can 
see from (0) that 5 runs quickly to a value such that the coefficient in Q changes its sign 
and ^13 increases again. Thus, the only case where ^13 can cross zero is the CP conserving 
case. This is interesting for future precision measurements of 6'i3, since the assumption 
of leptonic CP violation at any scale leads to the conclusion that the weak-scale value of 
6'i3 does not vanish. We illustrate the corresponding large effects in the evolution of 5 in 
Fig. m For all our plots, we use the software packages REAP and MPT associated with [33]. 

We can understand this feature also differently. In the above approximation, we can 
write = ^136""^ and by inserting the RGEs for 6*13 and 6, we find in the limit ^^13 

at 

- 2^"^;w'l"^ 2, ^in(2Msin(2M- (7) 

32 [m^ — ml) [m^ — 7712) 

For ^13 — i> we find thus that the RG change of Ues is along the real axis and Ue3 can 
therefore only become zero if it is real. The imaginary part of Eq. (0) allows to determine 
a minimal value of 6^13 as (6'i3)min — ^i3(/^) sin5(/i) where any scale /i can be used. 

Furthermore, let to denote the turning point of 6'i3 characterized by 6 = 7r/2. The 
'asymptotic' behavior 6(t — to) — —^ito — t) = vr — 6(to — t) is a consequence of the fact 
that 6 is an odd function in 6'i3^. This allows to understand why 6 approaches vr — S{mz) 
for large fi in Fig. [T] 

^The evolution of the weak basis CP invariant has already been studied in [25]. 

^This is in contrast to the case of Majorana neutrinos, where the memory to CP violation can be stored 
in the Majorana phases, and ^13 can cross zero even in the presence of leptonic CP violation [27,32]. 

^We have introduced tt in order to keep ^13 positive as we use the convention that ^13 is always positive, 
and a possible sign flip is absorbed in a jump of S. 
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13 J 
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Figure 1: Evolution of (a) sin^(20i3) and (b) 5 for small values of ^13. We use tan/3 = 50, 
612 = 0.1°, nil = 0.1 eV and best-fit values for all other parameters. The solid (black), dash- 
dotted (red) and dashed (blue) curve shows the evolution of (a) ^13 and (b) 5 for 6 = 45°, 20° 
and 1° at the electroweak scale, respectively. ^13 cannot become in any of the cases. 



The evolution of the mass eigenvalues is given by 



IGvr rhi 



cos 9i2 COS 623 sin 6'i3 + sin 612 sin 6 



723 



— - cos (5 sin6'i3 sin(2 6'i2) sin(2^23) 



+ ay>mi 



la) 



IGvr m2 



sin^ 6*12 cos^ 6*23 sin^ 6^13 + cos^ 6*12 sin^ 6*23 



+ -cos(5 sin6li3 sin(2^i2) sin(2 6'23) 
167r^ m3 = {C y1 cos^ O13 cos^ 6*23 + m-i . 



+ au} m2 



(8b) 
(8c) 



ay represents the flavor- independent part of the RGE, and is given in ()A.5d|l . Clearly, the 
dominant RG effect of the masses is a common rescaling governed by a^. In addition, for 
large tan/3 in the MSSM, there are corrections specific to the individual mj. In leading 
order in ^13, these effects tend to decrease /\m1^^ for a normal hierarchy and to increase 
Am^^jjj for an inverted hierarchy when running down. 

Fig. 121 shows an extreme example of the evolution of the mass eigenvalues and the 
corresponding Amg^i and Am^^^^. For large tan/5, there are substantial deviations from 
the flavor-independent scaling of the mass eigenvalues. The latter can be approximately 
inferred from the curve of mi in Fig. |21 (a). 
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(a) mi{iJ,)/mi{mz)- 



(b) Am? (/^)/Am? (mz). 



Figure 2: Example of the evolution of the mass eigenvalues and the Am^s. We choose nii{mz) = 
0.1 eV, 5{mz) = 50°, tan/? = 50 and a SUSY breaking scale of 200 GeV, and best-fit values 



3 Radiative corrections and future precision experi- 
ments 

An important question is if future experiments will reach a precision which allows to draw 
interesting conclusions from quantum corrections. There exist, for example, models where 
^13 vanishes at the GUT scale, but RG corrections still lead to a finite value of ^13 at low 
energies. A certain finite value of is therefore guaranteed unless the initial values at the 
GUT scale and the independent RG effects are adjusted to cancel each other. From the 
discussion of the previous section, we would also know that the CP phase 6 would vanish 
for all scales for Dirac neutrinos, while it could become finite for Majorana neutrinos. A 
finite value of 5 and 6*13 would thus exclude either Dirac neutrinos or 9i2,{Mgut) = 0. 
Similar arguments can be made for other quantities of the leptonic mixing matrix. ^23, for 
example, is within current experimental errors compatible with 45°. However, as in the 
case of 6*13, certain deviations are expected from RG effects even if 45° is exactly predicted 
at the GUT scale. Future precision measurements of neutrino oscillations may therefore 
allow interesting tests of flavor models and related renormalization group effects. 

Atmospheric neutrino data [34] and results from the K2K long-baseline accelerator 
experiment [34] currently determine ^m\^ = (2.2^04) x 10~^eV^ and 623 ~ 45° [34, 
35], whereas solar neutrino data [36-43], combined with the results from the KamLAND 
reactor experiment [44] lead to Am^^ = (8.2l°J) x 10"^ eV^ and tan^ ^12 = 0.39l[5:[J^ [35]. 
These results are to a good approximation still described by two independent two flavor 
oscillations. The key parameter for genuine three flavor effects is the mixing angle ^13 
which is so far only known to be small from the CHOOZ [45, 46] and Palo Verde [47] 



otherwise. 
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experiments. The current bound for ^13 depends on the value of the atmospheric mass 
squared difference and it gets weaker for Am|^ ^ 2 x 10~^eV^. However, in that region an 
additional constraint on 6*13 from global solar neutrino data becomes important [48,49]. At 
the current best fit value of Am|i = 2.2 x 10"^ eV^ the 3a bound is sin^ 613 < 0.041 [35]. 

Genuine three flavor oscillation effects occur only for a finite value of ^13 and estab- 
lishing a finite value of ^^13 is therefore one of the next milestones in neutrino physics. 
Leptonic CP violation is another three flavor effect which can only be tested if ^13 is fi- 
nite. There exists therefore a very strong motivation to establish a finite value of ^13 and 
then leptonic CP violation [50-54]. Different experimental projects are therefore under 
construction or are being planned in order to achieve these goals. It is useful to divide 
the future into what can be achieved with specific current or next generation projects and 
what may be achieved with long term projects. The MINOS [55] project, which started 
already data taking, and the CNGS projects ICARUS [56] and OPERA [57], which are 
completing construction can be considered as "current projects" . Beyond that exist other, 
more ambitious "next generation" projects like the superbeam experiments J-PARC to 
SuperKamiokande (T2K) [58] and the NuMI off-axis experiment NOi^A [59]. In addition 
there are "next generation" plans for new reactor neutrino experiments [60] with a near 
and far detector. A first interesting question concerns improvements of Am|]^ and sin^ 623. 
In Tab. Ewe show the relative precision which can be obtained in the future in comparison 
to the current precision, as obtained from a global fit to SuperKamiokande (SK) atmo- 
spheric and K2K long-baseline data [48,49]. We observe from these numbers, that the 
accuracy on Arrig^ can be improved by one order of magnitude, whereas the accuracy on 
sin^ 023 will be improved only by a factor two. 





|Am2j 


sin^ 623 


current 


88% 


79% 


MINOS+CNGS 


26% 


78% 


T2K 


12% 


46% 


NOz/A 


25% 


86% 


Combined 


9% 


42% 



Table 1: Relative precision of |Am|^| and sin^ ^23 at 3fT for the values Am|;^ = 2 x 10~^eV^, 
sin^ 623 = 0.5. The last row is the relative precision which can be obtained by combining all 
experiments (from [61]). 

Tab. n depends on the value of Am|]^ and the sensitivity suffers for all experiments 
for low values of Amg^^. T2K will provide a relatively precise determination of Am|^ for 
Amg^ > 2 X 10~^eV^. Although NOi/A can put a comparable lower bound on Amg^, 
the upper bound is significantly weaker, and similar to the bound from MINOS [61]. The 
reason for this is a strong correlation between Am|^ and 623, which disappears only for 
Arri^i ^ 3 X 10~^ eV^. From Tab. Q one can also see that only T2K is able to improve the 
current bound on sin^023- The main reason for the rather poor performance on sin^^23 
is the fact that these experiments are mostly sensitive to sin^26'23. This implies that for 
023 ~ vr/4 it is very hard to achieve a good accuracy on sin^ 6'23, although sin^ 26'23 can be 
measured with relatively high precision [62]. 

An important task of the next generation long baseline and reactor experiments of the 
coming years will be to establish the three flavored-ness of the oscillations. The sensitivity 
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to a finite value of the key parameter ^13 is only modest for MINOS, OPERA and ICARUS. 
Double Chooz, T2K and NOz/A can do much better. The sin^ 26'13-limits of the beam exper- 
iments are, however, strongly affected by parameter correlations and degeneracies, whereas 
new reactor experiments provide a "clean" measurement of sin^26'i3 [63]. Altogether these 
experiments will provide an improvement by about a factor ten for sin^ 26'i3 over the ex- 
isting limit. In addition, the KamLAND [44] (and solar neutrino) data will also further 
increase the accuracy for Am2i and Ou- An accuracy of 5% for Am^2 20% for sin^ 612 
is expected. Further improvements are possible, e.g. by loading the SuperKamiokande de- 
tector with Gadolinium, which might lead to an error of 3% for Am^i and 15% for sin^ 9i2, 
both at 99%CL [64]. 

Beyond the next generation accelerator and reactor based oscillation experiments exist 
much more ambitious projects like the JHF-HyperKamiokande project, beta beams and 
neutrino factories^. Such experiments clearly require further R&D before they can be 
built. However, assuming current knowledge, such facilities appear to be possible and they 
will lead to a precision at the level of percent or even below. With a neutrino factory, for 
example, a sensitivity to a finite value of sin^ 2^13 might be achievable down to 10~^. 

It is interesting to compare these perspectives with RG effects. To illustrate the RG 
effects, we start with initial values for the mixing parameters at the GUT scale, Mgut = 
3 X 10^^ GeV, assuming that these values find an explanation in a more fundamental 
theory.^ These initial values are then compared with the corresponding mixing parameters 
at TUz- In all our illustrations, we assume msusY = 1 TeV, and a normal mass hierarchy. 
The simple expressions (Eqs. (H}-© and (0)) allow a quick estimate of the RG effects. A 
more precise evaluation requires a numerical analysis for which we use the Mathematica 
package REAP [33], which is publicly available''. 

The mixing angles 612 and 623 turn out to be rather unstable for a degenerate spectrum 
(cf. Fig. 01). As a consequence, a Dirac version of quark-lepton complementarity [68-70] 
can only work for certain mass eigenvalues and ratios of the Higgs VEVs and tan/? (for 
the discussion of the RG effects in the see-saw Majorana case see [33,71,72]). This means 
stability of 612 is only given in models with hierarchical masses and/or small tan/?. Note 
also that for an inverted hierarchy 612 is unstable. This means that concerning 612 RG 
effects are in general an issue. RG corrections to the special value 623 = 45° can be 
comparable to the precision of upcoming experiments. Again, this happens for a quite 
degenerate spectrum and/or large tan/5. 

The running of ^13 depends crucially on its initial value. We illustrate this by plotting 
the radiative correction to sin^ 26*13 in Fig. |31 Most important is the second term in Eq. ^ 
which is dominant for not too small 613. As a consequence we find that, for ^13 = at 
the high scale, running in general does not generate a measurable value at the low scale. 
Only for the most optimistic sensitivities, a quite degenerate spectrum and large tan (3 
this conclusion can be avoided. On the other hand, if ^13 is not tiny, RG effects can be 
comparable to the precision of upcoming experiments (except for small tan/?). 

Finally, let us discuss corrections to 6. From the previous discussion in SecElit is clear 

^See [65] for a comparison and for references. 

®One could, for instance, enjoy the possibility of fixing the initial values by continuous (e.g. [66]) or 
discrete (e.g. [67]) symmetries. In this case, RG effects add to the corrections arising from the breakdown 
of those symmetries. 

''See [http : // www . ph . turn . de/ ~rge/| 
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(a) Arg6'i2 = 6*12(^2) - 6'i2(Mgut)- 



(b) Arg023 = |0.5 - sin2 023(MGUT)|. 



Figure 3: Radiative correction to (a) 612 and (b) 623 for 612 = 33°, ^13 = 10°, 623 = 45° and 
6 = 90° at = rriz as a function of tan/? and mi. These contours remain to a large extent 
unchanged when varying ^13 in the allowed range and 5 arbitrarily. 
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(a) ei3{mz) = 10°; S{mz) = 90°. 



(b) 013 (mz) = 0. 



Figure 4: Radiative correction to sin 29i3, defined as ArgsIu^ 26*13 = | sin^ 26'i3(Mgut) 
sin^ 20i3(m^)|, as a function of tan/3 and mi. We take 0i2{mz) = 33° and ^23("iz) = 45°. 



that small ^13 corresponds to an unstable configuration with large RG effects, even for a 
hierarchical spectrum (cf. Fig. ini(b)). In particular, RG effects are generically comparable 
with the precision of future experiments such as the combination T2K+NOi/A+Reactor-II, 
T2HK and NuFact-II (see [73] and references therein). 

Let us finally mention that RG effects for Dirac neutrinos will always result in a rescaling 
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tanjS tanyS 

(a) 9i3{mz) = 10° (b) 0i3(mz) = 0.1° 



Figure 5: Radiative correction to 6, defined as Arg<^ = \S{mz) — 5(Mgut)| for (a) diz{mz) = 10° 
and (b) di3{mz) = 0.1° as a function of tan/3 and mi. We use 5{mz) = 90°. 

of the mass eigenvalues. Beyond that, in the framework of the SM, mixing parameters are 
quite stable. The only exceptions are 9i2 for very degenerate masses, and 5 for tiny ^13. 
On the other hand, in the MSSM, RG effects are increased by tan^ /3, i.e. by up to three 
orders of magnitude. 



4 Summary 

Assuming Dirac neutrinos, we have derived renormalization group equations for leptonic 
mixing parameters. The results share several features with the corresponding equations for 
Majorana neutrinos. However, Dirac running is more predictive, as the Majorana phases 
are unphysical in this case. This makes it possible to specify the amount of renormal- 
ization group evolution unambiguously as soon as mi and 5 (and tan/?) are known. The 
renormalization group evolution alone does not yield an explanation of the largeness of the 
leptonic mixing angles (for an analogous and very clear discussion for Majorana neutrinos 
see [74]). Yet it may account for radiative enhancement of ^^12, and possibly also of ^23, 
since both can increase significantly in the MSSM when running down. 

Most importantly, we find that in phenomenological studies renormalization group 
effects for leptonic mixing angles can in general not be neglected. This can be understood 
from the fact that 6ij = f{mi,9ij,S)/{m'^ — m^) which becomes singular if rrii — > ruj and 
vanishes if the mixing angles are zero. We have thus traced back the relative enhancement 
of the quantum corrections of leptonic mixing parameters as compared to quark mixings 
to two reasons. First, the mass hierarchy which suppresses the renormalization group 
running, can be much weaker. Second, the mixing angles are larger so that the parameters 
are further apart from the renormalization group stable situation where all mixings are 
zero. 
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As there is no suppression of the running by phases, the renormahzation group correc- 
tions should in general be taken into account even for a strong hierarchy to accommodate 
the precision of future experiments. Renormahzation group corrections are especially rel- 
evant if the mass spectrum is non-hierarchical, and tan/3 is large in the MSSM. Hence, 
similarly to the case of Majorana neutrinos [27], also in the Dirac case the non-observation 
of deviations of the angles from special points, e.g. of 612 from n/A — 1^12 (with 1^12 being 
the Cabibbo angle), of ^13 from and 623 from 7r/4, may restrict the parameters such as 
the absolute neutrino mass scale. The current experimental data already has the necessary 
precision to indicate disfavored parameter regions. It may also point to exactly realized 
symmetries and our formulae can hence be used to identify possible symmetries. Whenever 
a symmetry is exact and fixes some mixing parameters, those mixing parameters have to be 
stable under quantum corrections. For instance, it has recently been pointed out [75] that 
for Majorana neutrinos and an inverted hierarchy the configuration ms = ^13 = is stable. 
Prom the analytic expressions it is obvious that this statement also applies to the Dirac 
case. Likewise, a quick inspection of the RGEs excludes most of the proposed symmetries 
from being exact. Our formulae are basis-independent, and thus allow to understand cer- 
tain features of the underlying theory, such as symmetries, in a basis-independent way. We 
have discussed this only for the case of CP symmetry, but it is obvious how the analysis 
can be carried over to other symmetries. In this context, it would be interesting to see 
if infrared fixed points with large mixings, as discussed in [74,76], can be obtained for 
(non-standard) Dirac neutrinos as well [77]. In this case, one may hope for a scenario 
where the large mixings are a consequence of running, and the mechanism of generation 
of neutrino masses is still related to the scale where gauge couplings meet. 

We conclude that in the light of future precision experiments, flavor physics might en- 
ter into an era of "precision model building". It seems possible to determine the mixing 
parameters to a remarkable accuracy, precise enough such that flavor models and the cor- 
responding renormahzation group effects become to a certain degree distinguishable. For 
a specific parameter and a desired accuracy, our formulae allow to estimate the renormal- 
ization group effects, and to judge to which extent a numerical analysis is required. 
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A Mixing parameters RGEs for Dirac masses 



A.l Lagrangian 



The Yukawa couplings are given by 



in the SM extended by right-handed neutrinos where 6 = ia2(if ■ In the extended MSSM, 
the Yukawa couphngs are analogously defined in the superpotential by 



The left-handed lepton and quark doublets are denoted by ^L and Ql, respectively. We 
assume that there is no Majorana mass term for the right-handed neutrinos. 

A. 2 /3-functions 

The relevant /^-functions for the down-type quark, up-type quark, charged lepton and 
neutrino Yukawa coupling matrices Y^i Y^, Y^ and Y^, read at one-loop [78,79] 



(A.2) 



(AnfYa = Y^ 



{c',YjY, + C^Y^Y^ + aa} , 



(A.3a) 




(A.3b) 

(A.3c) 
(A.3d) 



(47r)2n 
(47r)2i; 



Y,{C:Y^^Y, + C:YX + c^i} , 
Y, {CtYX + C:YX + o^A , 



where 




(SM) 
(MSSM) 



(SM) 
(MSSM) 



(SM) 
(MSSM) 



(SM) 
(MSSM) 




(SM) 
(MSSM) 



(SM) 
(MSSM) 



(SM) 
(MSSM) 



(SM) 
(MSSM) 



(A.4d) 



(A.4b) 



(A.4a) 



(A.4c) 
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and 



-bl-lgl-^gl + TsM. (SM) 
~ ^ 3iv(yiy,) + iv(y,ty,)-^^f-3^i-f ^1, (mssm) "^^^ 



-M-lal-^gl + TsM. (SM) 

Tr(FjF,) + 3 Tr(Fji;) - If ^2 _ 3 ^2 _ ^ gl , (MSSM) 

-!^?-!^i + TsM, (SM) 



(A.5b) 



^ 3iv(yJye) + TV(yJi;)-|^?-35i , (mssm) 

= . -i^n^i+rsM, (SM) 

Tr(yJn) + 3Tr(yJyj-|^f-3^|, (MSSM). 

Here, we define TgM = Tr y^ty^ + yjy^ + Sy^y^ + 3yjy„ , and use GUT normafization 
for gi. 

A. 3 General derivation 

In this subsection, we will perform a general analysis applicable for any Dirac masses, and 
treat the evolution of lepton and quark masses and mixings only as a special case. 
We derive the running of mixing parameters for a RGE of the form 

IGtt^-^// = F^H + HF + fH, (A.6) 

vJ.6 

where / is real and H is Hermitean, so that we can diagonalize it in a 'reference basis', 

-H-U = D . (A.7) 

In the application in the main part, F corresponds either to C^yjy^ + C^Y^Y^ (or 
C^YjYe + Ci^YjY^ for the lepton sector), and H to Y^Ya (or YjY^). The 'reference basis' is 
the basis where Y^Yu (or Y^Ye) is diagonal a,t t — to. U denotes then to the CKM matrix 
f^CKM (or the MNS matrix C/mns)- / denotes the diagonal parts of the /3-function, / = 2ad 
(or / = 2a^). 

Now we perform an analysis very similar to what is done in [27] which is based on 
[32, 79, 80]. We can differentiate the relation H -D-U^ 

^{U-D-U^) = if ■ D -U^ + U ■ D ■if'' + U ■ D ■U'' 

= -^(F^-U-D-U^ + U-D-U^-F + fU-D-U^) . (A.8) 

Multiplying by W from the left and by U from the right yields 

U^-U-D + D-U^-U + i) = -^(F'^ -D + D-F' + fD) , (A.9) 

where F' ^U^^ -F-U. For the quark case, F' ^CjD + C^ U^Y^Y^ U. We will see below 
that only the off-diagonal components are relevant for the RGEs of the mixing parameters. 
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The evolution of U can be written as 



^f/ = U-X , (A.IO) 
where X is anti-Hermitean. Inserting this relation yields 

D + X-D + D-X^ = -^(F'^ -D + D-F' + fD) , (A.ll) 
or, by using the anti-Hermitecity of X, 

b = if D + F'^ ■ D + D ■ F') - X ■ D + D ■ X . (A.12) 

Denoting the entries of D by yf, i.e. D = diag{yf, yl, ?/|), we find 

j^yl = Yl^[fy^ + iF:: + F:^y^] , (a.is) 

i.e. the terms proportional to X have dropped out. This equation corresponds to a RGE 
for the running mass eigenvalues, defined by mi{t) = \yi{t) \ v with v fixed, of the form 

{ATifmi = iReF^- + a)mi. (A.14) 

By analyzing the off-diagonal parts we obtain 

X,, - X,, y] = - ^ [(F't),^. y] + y] F[^ . (A.15) 

This can be converted into equations for real and imaginary part of X, which, since F is 
Hermitean, can be combined to 

levr^X,, = 4^ 4 • (A.16) 

The diagonal parts of X remain undetermined. However, this is not a problem, since 
they only influence the RG evolution of the unphysical phases.^ 

So far, we have derived the differential change of the CKM matrix due to the RG 
corrections for Y\Yd (cf. Eq. ()A.10|) ). In the Majorana neutrino case, the analogous differ- 
ential equation already describes the evolution of the MNS matrix since Y^^Ye doesn't get 
rotated by the RGE.^ For Dirac neutrinos, Y^Yf, gets rotated only by terms proportional 
to the squares of Dirac Yukawa couplings, hence those rotations can safely be neglected. 
In the quark sector, the radiative rotation of Y^Y^ represents an important effect, as we 
will argue in the following. 

^Note that the Majorana phases are unphysical in the the Dirac case as well. 

^This is only true at leading order in where M denotes the scale of the effective neutrino mass 

operator (e.g. the see-saw scale) [81]. 
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A. 4 Contribution from the change of Yu 

Here, we specialize to the quark sector as the analogous effect is irrelevant for Dirac 
neutrinos. 

The RGE for contains non-diagonal terms so that continuous re-diagonalization is 
required. Since the mixing matrix Uckm is defined as the matrix which diagonalizes yJYii 
in the basis in which Yu is diagonal, Uckm receives an additional contribution from the 
running of Yu, 

= ?7cKM ■ -'^ + term stemming from the change of . (A. 17) 

dt 

To evaluate this change, we can essentially repeat the steps of the previous subsection. 
Introducing a matrix U which diagonalizes Y^Yu in the reference basis (implying U{t = 
to) = 1), i.e. 

U^Y^U = diag{yl,yl,yl) , (A.18) 
we arrive at 

= U-X, (A.19) 

where the off-diagonal entries of X are given by 

WX,, = I^F,,. (A.20) 
ilj ill 



Completely analogous to IA.31 

F' = -F-U , (A.21) 
and at t = to 

F' = Cf^U DU^ + ClY^Yu. (A.22) 
Again, only the off-diagonal terms influence the RGEs of the mixing angles. 

A. 5 Mixing parameter RGEs in the quark sector 

As f/cKM = U-^U = UW, the RGE for the CKM matrix reads 

^f^CKM = X^ ■ Uckm + Uckm ■ X . (A. 23) 
To proceed, we label the mixing parameters by 

{6} = {dl2, Oi3, O23, 5, 4, ^r, <^i, ^2} , (A.24) 
and evaluate the derivative of Uckm, 

Uckm = f/cKM fe}, fc}) • (A.25) 
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Observe that the resulting expression is hnear in ^fc. By solving a system of linear equations 
of the form 



J2 4x ik + i 4x ik = Rx, (A.26) 



k 

where 



Rx = Umns ■ T + ■ Umns , (A.27) 

we thus obtain a set of linear equations for the C,k- RGEs for the matrix elements have 
been derived in refs. [79,80]. 

From these, we obtain the RGEs for the mixing angles in the quark sector. Neglecting 
all Yukawa coupling except for yt and yb, they read 

^12 = %^ cos(i9i2) { [(3 - cos 2 T^is) cos 2^23-2 cos' i9i3] sin l9i2 

+ 4 cos^cp cos'(9i2 sin-i^ia sin2'(923} , (A. 28a) 

^13 = "g^y^' [2 Ct yl + Q yl{l + cos 2 ^23)] , (A.28b) 

^23 = "g^y^' K (1 + cos 2 ^13) + 2 y,'] . (A.28C) 

It turns out that finite ys and yc corrections yield an important but sub-dominant effect 
for di2. The dominant term in the RGE of 5cp is 

(ju 2 2 

'^cp = ^Vr^^Hv cos^9i2 cos7?23 sin5 sin79i2 sin^923 X ^9^3^ . (A.29) 
8 vr2 {yl - y2) 

A. 6 Mixing parameter RGEs in the (Dirac) neutrino sector 

In order to derive analogous RGEs for the leptonic mixing parameters, observe that the 
RG change of Y^Ye is quadratic in neutrino Yukawa couplings, i.e. strongly suppressed. 
Thus, we can safely neglect the X contribution, 

-rUyws = ■ f^MNS + f^MNS ■ AT ~ Umns ■ X , (A. 30) 

dt 

where X is now related to F' by Eq. (|Xl6|l . and F' = Q Li + C^ f^MNS^e"^^e t^MNS at t = to- 
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